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Abstract 

' In order to investigate the evolutionary process of many deterministic dynamical systems with unfixed 

TJ*± , parameter, a set of dynamical models with parameter changing continuously and the accumulation of 

this change might be large is introduced and discussed. The boundary crises and the period-doubling 
bifurcations are found suppressed and scaling properties for these phenomena are exploited. Due to this 
suppression, the period-doubling bifurcations seem no longer continuous. We further discuss the possible 
■ applications of these models. 

> : 

"nF ■ PACS number: 05.45. +b 

o ; 

Introduction. One of the practical methods to investigate the phenomena common in many scientific 
fields such as physics, biology, ecology, geophysics and economy, etc is by modeling them either with dif- 
ferential, difference-differential, or pure difference equations. After the models have been studied, these 
Q\ [ phenomena can be understood and even be predicted. In this Letter, we will introduce a model which can be 

used to model and then study the behavior of many natural systems with unfixed environment parameters. 

The motivation of introducing such models is our recent carefully examination on the work which was 
tempted to apply the chaotic dynamics to understand and predict the natural phenomena [see, e.g. ref 1-5]. 
C. Nicolis and G. Nicolis proposed some evidence for the climate which caused a debate on whether climatic 
attractors exist or not 1 . The correlation dimensions 6 for many economical data sequences had also been 
calculated which were taken as the evidence that the economical chaos might exist [see, e.g., ref. 2]. And 
recently, the short term prediction algorithm was used to test the predictability of some economical data 3 . 
Chaos in heart and brain was also reported and tempted to be controlled 4 . All of these studies are based on 
the assumption that these natural systems are environment independence, or the data of these dynamical 
systems are take in a limited time interval so that accumulation of the environment parameter change is 
sufficiently small. Then the familiar concepts such as phase-space, bifurcation diagrams, orbits, attractors, 
return maps, basins of attractor, manifolds, etc. can be used to extract information of these natural systems 
7 ~ 8 . In fact, many of the natural systems such as all biological, ecological and economical systems have 
developed into other ones before they reach any asymptotic states. The environment parameters change so 
quickly that accumulation of the parameter change is significantly large in a limited time interval which 
has changed the dynamical behavior of the systems. The data sequences from these systems, on one hand, 
are some kinds of transients in short term interval (so that the accumulation of the parameter change is 
small), and on the other hand, are the sets of points in different environment parameters in a relative long 
time interval. It is clear that the above listed techniques are applicable only when the change of environment 
is negligible. Thus, this assumption not only restricts the study on the global behavior of these natural 
systems, but also might cause misunderstanding. 

In order to study the dynamical behavior of these natural systems, and clarify the possibility and limi- 
tations to apply the familiar concept to study natural systems, in this Letter we will study the dynamical 
models with parameter changing continuously and the accumulation of this change might be large. 
Explicitly, we will study the models in the form 

X n+1 =F(X n ,P n ), (1) 

P n+1 = G(P n ) + P n . (2) 
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where X n and P n represent one-, two- or high-dimensional variable and parameter vectors, respectively. 
Considering the continuous evolutionary process of the natural systems, G(P n ) is usually small (If Eq. (1) 
represents a difference form of differential equations with very small time step St, G(P n ) — > as St — > 0). 
One particular case for this model has already been discussed in detail which is the periodic forced driven 
systems. Some of the special cases of this model with the parameter close to the bifurcation points or G(P n ) 
is a noise had also received attention recently 9 ~ n . In the present Letter, we will concentrate on the situation 
with the parameters change unperiodically and ever irregularly, and the accumulation of the change of the 
parameters might be large (hereafter wc call these systems the structure-variable systems, namely, SVS's) . 
We will discuss the behavior and possible applications for these models. 

Due to the unperiodic change of the parameters, the dynamical systems will never approach asymptotic 
states. All the data we observed are transients. When the change of the parameters is not too fast (G(P n ) 
is sufficiently small), for a short observation time interval, the SVS represented by eqs. (1-2) might exhibit 
some behavior similar to those of Eq. 1 with fixed parameters P n . For P n close to the parameters for 
bifurcations or boundary crises 12 of Eq. 1, the behavior of the SVS (1-2) is significantly different from 
that of system (1) with fixed parameters. The crises and period-doubling bifurcations are suppressed and 
scaling properties about these phenomena can be found. Due to this suppression the periodic orbits and 
period-doubling bifurcations are difficult to be observed directly. 

Example For a detail idea of the behavior of these SVS's we take the following map 

J %n+l /(Mto; ^n) — 1 l^n^ni (o\ 
\ Mn+1 = Hn + Sfl, 

as an example. Fig. 1 shows a typical case for a trajectory of this map. This diagram shares some 
similarity with the bifurcation diagram of the Logistic map in the same parameter interval. Unlike that in 
the bifurcation diagram, no point is thrown away as transients. In fact, this diagram reflects the evolutionary 
process of the dynamical systems described by the Logistic map (we will call it the developing diagram 
hereafter), which corresponds to an idea model of the insects population on an island 13 with the environment 
of the island changes (The area might shrink because of the housegreen or other effects) . The value of the 
parameter step Sfi reflects the rapidity of evolutionary process of the dynamical system. It is clear that 
when the width of any periodic window is small or comparable with the value of 8/j,, this periodic window 
disappears in the developing diagram. 

A direct observation of the Fig. 1 is that the period-doubling bifurcation seems no longer a continuous 
one. Comparing to the bifurcation diagram, the parameter of period-doubling is shifted. Now we estimate 
the rapidity of the convergence of the period-doubling bifurcations and the shift with respect to Sfi in these 
developing diagram. 

Consider a period-doubling bifurcation occurs at fi c in the bifurcation diagram and at in the 

developing diagram with a parameter step <5/i. /i is the controlling parameter. ^id(Sn) > In fact, from 
the parameter fi — fi c , the developing diagram goes alone the lose-stability periodic orbit until it reaches 
fi = /id, due to the slowing down of the convergence and the inertia. 

Consider a period p orbit which should be a certain fixed point x* of the following map 

Xn+p = / O / O • • • O f(fj,, X n ) . (4) 
v v ' 

p times 

Let 

\e n \ diminishes as e -71 !'* - ^! for each iterate of the map (4) at parameter fi, and A is a positive constant 14 . 
Consider that we go forward k steps from fj, = /i c of map (4) and reach jia = /i c + kSfj,. Assume that the 
initial deviation of the period p orbit is e^°\ and the final deviation is e^ k \ 

| e (fe)| a e - A ELi lw-^h e (o)i = e-^ELiWleC)! = e -^(k+i)l«Ml/2| e (o)i 

For all small \5fj,\ we assume |e( fc )/e(°)| to be an approximate constant, so that k(k + l)Sfi is approximately 
a constant. Reminding that kSfi = (id — fi c , we obtain that 

\Hd(Sfi) -Mc| oc y/\Sfl\. (5) 
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We have checked this relation numerically for many period-doubling bifurcations in many models. This 
relation holds approximately provided the parameter step S/i is not too large. 
At \x = Ha the rapidity of convergence is proportional to 

The larger the value of \jid — A* c j > the smaller the value of e _j4 l Al<i_ '* c l. Thus the larger the value of 8/i, the 
more rapid the convergence (see Eq. 4). This makes the period-doubling bifurcation no longer a continuous 
bifurcation for not too small 5fi. 

The parameter for a boundary crisis is also shifted as shown in Fig. 2. In a similar way we get the scaling 
law for the parameter shift A/i of the crisis with respect to the controlling parameter step Sn as following 

A/iocd/i/X, (6) 

where A is the maximal Lyapunov exponent. 

The parameter shifts for both the boundary crisis and period-doubling bifurcations cause the so-called 
"artificial" hysteresis that the trajectory is still stay in a lose-stability attractor while another attractor has 
already existed. 

Besides the enviroment parameter changing of of many dynamical systems, it should also pay attention 
on the method of extracting data from the natural systems, in which we usually use the average or sliding 
average values such as the economical data per hour, day month and even year, and the 6-values in seismic 
events which are the sliding averages of seismic events in a definite time interval 15,16 . The periodic orbits 
from these data are totally obscured. In Fig. 3 we shown the sliding average values of the SVS (2). It is 
clear that the periodic orbits correspond to some dips. This gives us a hint that a tangent bifurcation might 
happen at this parameter range and periodic orbits might exist. 

Possible applications There are various applications of the developing diagram. Recently, short term 
predictability 17 on deterministic dynamical systems with fixed parameters has been obtained a great deal 
attention. Are the data from a dynamical system with changing parameter predictable? With our model, 
we can show that that the short term predictability is possible for these natural phenomena no matter 
the environment parameter changes slowly or quickly, provided that the changing of the parameter is 
governed by some rules (Though we might not know exactly what they are, see detail in our recent work, ref, 
18). With this idea, we can understand that it is reasonable that some economical data 4 and the b- values 19 
for seismic events are predictable. 

The developing diagram can also be used to clarify the real physical meaning of the various character- 
izing quantities (correlation dimension, Lyapunov exponent, etc) which were calculated for many natural 
dynamical systems 18 . In fact, due to the changing of parameters, the reconstructed figures for the data from 
a developing diagram is very similar to that for a noised dynamical system. 

An important application of the developing diagram lies in the study of the subsystems (or open sys- 
tems) of dynamical systems. The influence of the other system can be considered as one or more changing 
parameters, which might be written by some equations phenomenally. This can help us to understand the 
restrictions of applying various perturbation and truncation theory to nonlinear systems since a small chang- 
ing of the environment parameter (by adding high order perturbations or truncations) will suppress or induce 
chaos, crises 20 . 

Summary and discussion In summary, we have introduced a set of dynamical models with parameter 
changing continuously and the accumulation of this change might be large. The parameters for boundary 
crises and period-doubling bifurcations are shifted and scaling properties for these phenomena are exploited. 
Due to these shifts, the period-doubling bifurcations seem no longer continuous. We also discuss other 
possible applications of the models. 

There are still a number of interesting questions remained to be investigated. Since the data series 
from our model (as well as many natural systems) are the set of transients for different parameters, how 
to extend the general mcthods[see, eg. 21,22] to extract information from the data series and then control 
and predict them is needed. The "artificial hysteresis" clearly has closely connection with many biological 
phenomena. The average values, also introduced in this Letter, should be paid more attention even for that 
from a dynamical systems with fixed parameters. We are afraid that a great part of noise we considered in 
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many natural systems results from the averaging process and the changing parameters. These questions and 
the other applications of the developing diagrams are now undertaken by the author and his collaborations. 
Some of them will be presented in an extended paper. 
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Figure Captions 

Fig. 1 The developing diagram of the logistic map in the parameter interval fi E [1.6, 1.82] with a parameter 
step Sfi = 10~ 5 . The initial point xq at fi = 1.6 is on the attractor. 

Fig. 2 Same as Fig. 1 with fj, e [1.999, 2.003]. 

Fig. 3 The sliding average values of the Fig. 1 with 500 data a moving. 
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